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Active flow control technology is finding increasing use in aerospace applications to control flow separation 
and improve aerodynamic performance. In this paper we examine the characteristics of a class of fluidic actua- 
tors that are being considered for active flow control applications for a variety of practical problems. Based on 
recent experimental work, such actuators have been found to be more efficient for controlling flow separation in 
terms of mass flow requirements compared to constant blowing and suction or even synthetic jet actuators. The 
fluidic actuators produce spanwise oscillating jets, and therefore are also known as sweeping jets. The frequency 
and spanwise sweeping extent depend on the geometric parameters and mass flow rate entering the actuators 
through the inlet section. The flow physics associated with these actuators is quite complex and not fully un- 
derstood at this time. The unsteady flow generated by such actuators is simulated using the lattice Boltzmann 
based solver PowerFLOW®. Computed mean and standard deviation of velocity profiles generated by a fam- 
ily of fluidic actuators in quiescent air are compared with experimental data. Simulated results replicate the 
experimentally observed trends with parametric variation of geometry and inflow conditions. 


Nomenclature 


D 

divergence angle, degrees 

lb 

pound 

m 

meter 

mdot 

mass flow rate, lb/sec 

mm 

millimeter 

k 

turbulence kinetic energy 

psi 

pound per square inch 

s, sec 

second 

t 

time, seconds 

V 

velocity magnitude 

Vavg 

time-averaged velocity magnitude 

V'rms 

root-mean square of perturbation velocity 

Vtotal 

instantaneous total velocity 

x,y,z 

Cartesian coordinates 

3-D 

three-dimensional 

e 

dissipation length scale 
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Abbreviations: 


AFC 

CFD 

DNS 

ERA 

LBM 

VR 


active flow control 
computational fluid dynamics 
direct numerical simulation 
environmentally responsible aviation 
lattice Boltzmann method 
variable resolution 


Superscript: 

o angle in degrees 

' perturbation quantity (e.g. V' = V-Vavg) 


Subscript: 

0 initial value 


I. Introduction 

There has been a growing emphasis on reducing the fuel bum of commercial aircraft to meet the stringent restrictions 
on emissions and noise imposed by regulatory agencies. A significant effort is being devoted to such activities at NASA 
through the Environmentally Responsible Aviation (ERA) program. One area of research under this program is to 
increase the lift/drag ratio and/or reduce the size of control surfaces on aircraft by making use of active flow control 
(AFC) devices, which will decrease the overall weight of the aircraft. The AFC devices have to be designed to maximize 
performance with the least amount of weight penalty and power requirements. In addition, these devices must have 
acceptable fault tolerance. 

Recent advances in computational fluid dynamics (CFD) methods have made it possible to simulate the effect of 
active flow control devices on major aircraft components, as evidenced by a series of publications in the literature. 1-5 
The paper of Collis et al. 4 summarizes the various issues encountered in active flow control applications from both 
theoretical and experimental points of view. Until recently, flow control was achieved primarily with synthetic or pulsed 
jet actuators. All of these devices consist of moving parts and controllers, which can cause operational problems in 
flight. 

During the last several years, fluidic actuators have come to prominence as viable flow control devices for improving 
the aerodynamic performance of high lift configurations by mitigating separated flow regions. Such devices do not 
contain any moving parts and rely primarily on the feedback mechanism triggered by the unsteady flow inside the 
actuators. DeSalvo et al. 6 and Woszidlo and Wygnanski 7 have examined the performance of such actuators in recent 
papers. Based on the available research, it appears that the fluidic actuators that produce span wise oscillating or sweeping 
jets offer a reliable and efficient mechanism for active flow control applications. 

CFD simulations are used here to examine the internal and external flow structures created by fluidic actuators 
operating in quiescent air. Numerical studies are conducted to quantify the effect of geometric shape changes and input 
conditions on the performance characteristics of such actuators. The computed solutions are compared with available 
experimental data for several configurations. 


II. Configurations 

The fluidic actuators under consideration here are closely related to the ones examined by Gokoglu et al. 8-10 The 
primary difference in the current configuration is the absence of the divider section used in their work to create twin 
outlets for the fluid to exit. In the current configuration, which is simpler to fabricate and is being used in several 
practical applications, the jet comes out of a single outlet, as sketched in Fig. 1. Here the x-coordinate aligns with the 
primary jet flow direction, whereas the z-coordinate refers to the spanwise direction. The y-coordinate is out of the paper 
in this figure and is normal to the x-z plane. Some of the geometric details are suppressed here (covered by gray regions 
in Fig. 1) because of the proprietary nature of the device. 

The inflow boundary is connected to a high pressure source to inject a constant supply of air into the actuator. The 
pressure differential between this air supply and the ambient atmosphere sets up a constant mass flow through the inflow 
section. The mass flow going through the actuator is directly proportional to the pressure differential, or gauge pressure. 
The flow accelerates up to the first nozzle before entering the main chamber, which contains constrictions placed on 
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the left and right sides, sub-dividing the main chamber. Due to the shape of this geometry, an oscillating unsteady flow 
develops inside the chamber such that in addition to the flow exiting through the second nozzle near the exit plane, 
there exists a backflow moving between the left and right side of the chamber. As a result, the flow exiting this nozzle 
oscillates from left to right in a cyclic pattern. The frequency and the span wise extent of the oscillatory jet depend on 
the precise shape of the actuator and pressure differential between the inflow boundary and the ambient atmosphere. 

We consider two types of fluidic actuators in this paper. The first type are formed with straight edges, sharp comers, 
and planar surfaces. The divergence angle of the exit nozzle for the baseline configuration, shown schematically in 
Fig. 1, is 124°. The divergence angle of exit nozzle for the second actuator of this type under consideration is 90°. Flow 
characteristics of similar actuators have been examined by Gokoglu et al. in a series of papers. 8-10 These will be referred 
to as Type I, 124D and 90D actuators, or for simplicity as 124D and 90D actuators. 

We also examine another type of fluidic actuator in this paper where the internal chamber consists of smooth curved 
edges instead of sharp corners and edges, as shown in Fig. 2. Woszidlo and Wygnanski discussed the geometric details 
and origins of such actuators in a recent paper. 7 At the start of the process ( t = to X the jet attaches to either wall of 
interior cavity, shown schematically in Fig. 2. The curving of the jet increases the pressure at the feedback channel inlet 
and induces reverse flow in the feedback channel that pushes the jet to the opposite side at a later time, t = to + dt. The 
process repeats itself at a frequency set by the precise shape of the actuator and the inlet conditions. A series of such 
actuators were embedded in a high-lift configuration consisting of a NACA 0021 airfoil with variable flap geometry. 
The effectiveness of these actuators to control the separated region for such airfoils was clearly demonstrated in Ref. 7, 
where significant improvement in measured lift was achieved with flow control compared to the baseline configuration. 
These actuators will be referred to as Type II curved- wall actuators, or as just curved- wall actuators for simplicity. 

III. Experimental Setup 

The fluidic actuators described in the previous section have been manufactured at NASA Langley Research Center 
and are being tested in a bench-top environment. A photograph of the curved- wall actuator model under consideration is 
shown in Fig. 3. A high pressure air supply is attached to the inlet boundary and the flow exits from the outlet face into 
quiescent air of the room. The mass flow entering the actuator and the magnitude of the exit velocity depend strongly on 
the pressure differential between the inflow plane and the ambient atmosphere. The mass flow rate entering the actuator 
is measured with a flow meter placed inline with the air supply tubing. A static pressure probe placed inside the inlet 
section is used to measure the inflow pressure. Hot-wire measurements are taken at several locations downstream of the 
exit plane. The raw data is processed to generate a time-series of the velocity signal. In addition to the instantaneous and 
mean velocities, the root mean square (rms) or standard deviation of the perturbation velocities are also extracted from 
the hot-wire data. 

Fluidic actuators of different depths have been fabricated and tested at NASA Langley. However, the results presented 
in this paper correspond to the actuators that are 0.25 inches (6.35 mm) deep. The throat width of the exit nozzle is nearly 
equal in size to the depth for these actuators. Therefore, the aspect ratio of the rectangular throats for these actuators is 
approximately one. 


IV. Simulation Method 

The numerical simulations presented in this paper were performed using the commercial software PowerFLOW® 
(version 4.3) which is based on the three dimensional 19 state (D3Q19) lattice Boltzmann model. The lattice Boltzmann 
model (LBM) can be derived from the velocity space discretization of the general Boltzmann equation, and it recovers 
the full unsteady Navier-Stokes equations. 11-14 The lattice Boltzmann equation can be solved in a DNS mode, where all 
of the turbulent scales are spatially and temporally resolved. However, for most engineering problems at high Reynolds 
numbers, the simulations are usually conducted in conjunction with a hybrid turbulence modeling approach where the 
small scales are modeled using a modified two-equations k — e turbulence model and larger energy containing scales 
are directly resolved. The turbulence model uses a swirl correction to locally reduce the influence of the modeled 
eddy viscosity to allow vortical structures to develop and persist without added artificial damping. 15-17 In addition, a 
wall model relying on algebraic general relations such as the logarithmic law of the wall, and extended to account for 
pressure gradients effects is used to capture the boundary layer behavior and relax the grid spacing requirement near 
the solid walls. 18 The turbulence modeling approach together with the inherently unsteady nature of the LBM ensures 
accurate and efficient simulation of the large structures especially in regions of separated flows and shows conceptual 
similarities with the recent DDES approach of Spalart. 19 This was the approach used for obtaining numerical solutions 
presented in this paper. 
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The lattice Boltzmann equation is solved on embedded Cartesian meshes, which are generated automatically within 
the flow solver for the configuration under consideration, irrespective of the geometric complexity. Furthermore, vari- 
able refinement regions (VR) can be defined to allow for local mesh refinement of the grid size by successive factors 
of two in each direction. 20 The PowerFLOW® code scales well on highly parallel clusters using up to thousands of 
processors, making it ideally suitable for large scale applications. The LBM used here has been extensively validated 
for a wide variety of applications ranging from academic cases using DNS 21 to industrial flow problems in the fields of 
aerodynamics, thermal management, and aeroacoustics, see e.g. Refs. (5, 17,22,23). 

A cross section of the computational grid (lattices) for the Type II, curved- wall actuator used in the PowerFLOW® sim- 
ulations is shown in Fig. 4. The 3-D computational domain consists of hexahedral voxels formed by extruding the 2-D 
grid in the y-direction. The volumetric mesh is very coarse in the outer region, which is essentially a quiescent (no flow) 
region and becomes successively finer as the computational boxes get closer to the actuator. The finest grids are used 
near the solid walls to resolve the high gradient viscous regions. The grid for this case consisted of a total of 29.7 million 
voxels (cells) distributed among 6 VR regions. Similar grids were created for the Type I actuator simulations. 

V. Results 

The unsteady 3-D flow fields created by the fluidic actuators exhausting into quiescent air were computed using the 
The PowerFLOW® code. The baseline conditions chosen for this work correspond to an inlet mass flow of 0.015 lb/sec 
at room temperature. Viscous wall conditions were imposed on all of the solid walls that form the actuator geometry and 
the internal components. A back pressure of 14.7 psi was imposed at the outflow boundary. Frictionless wall boundary 
conditions were imposed at rest of the computational boundaries. 

The computations were run for large number of time steps to generate a database from which statistically converged 
mean and perturbation flow quantities could be extracted. Instantaneous flow field data was collected at various locations 
corresponding to experimental measurement locations downstream of the actuator exit. Computed time-averaged mean 
and perturbation velocities for three different actuators are compared with the measured data. These results are discussed 
in the following sub- sections. 

A. Type 1, 124D actuator 

The first test case chosen for numerical simulations is the baseline 124D actuator which produces a sweeping jet that 
oscillates in the spanwise direction, z. To compare the performance and characteristics of such actuators with the ones 
that produce planar (non-sweeping) jets, the baseline configuration was modified as follows. Two 10 mm long solid 
blocks shown schematically as blue rectangles in Fig. 1, were inserted between the external walls and the internal diverter 
blocks on both sides to block the feedback channel. Although the flow field associated with this modified configuration 
is still unsteady in nature, the insertion of these solid blocks had the desired effect of creating an exit jet that resides 
mostly near the mid-section, instead of oscillating in the spanwise direction. The velocity magnitude contours from the 
numerical simulations for the modified configuration are shown in Fig. 5. Notice that the jet stays close to the center 
for this case. These additional solid blocks were then removed to recover the original baseline actuator configuration. 
The computational results for the original 124D configuration verify the presence of sweeping jet oscillating from left 
to right. Instantaneous velocity contours at two different times from these simulations are shown in Fig. 6. Results in 
Fig. 6 (a) indicate that most of the flow exiting the nozzle of the actuator at this instant in time is leaning towards the left 
side, thus leaving the flow on right side essentially undisturbed. The situation is exactly reversed in Fig. 6 (b), where the 
jet is leaning to the right. Time animation of these results (not shown here) is even more helpful in demonstrating the 
underlying flow patterns and the highly unsteady nature of the flow for these cases. 

The time-averaged and root mean square (rms) perturbation velocity profiles at the x=6 mm location for these cases 
are compared with the experimental data in Fig 7. The perturbation velocities shown in this figure have been normalized 
for each data set with the corresponding maximum average velocity. For the configuration with blocked feedback chan- 
nel, the maximum value of time-averaged velocity is almost 50% larger compared to the baseline 124D configuration. 
The opposite trend is observed for the normalized perturbation velocities (V'rms) in Fig. 7 (b), where the peak values 
for the blocked channel configuration are seen to be nearly half of the values attained for the baseline configuration. 
The predicted velocity peak value for the original 124D configuration is slightly higher compared to the measurements, 
while the situation is reversed for the blocked channel case. Nonetheless, the overall agreement of the computed average 
velocity profiles with measurements is quite good. Note that the predicted velocity profiles are fuller (broader) in shape 
for both of these cases compared to the experimentally measured profiles. 

Although the perturbation velocities are under-predicted for both cases, the computed profiles replicate the shape and 
the most salient features of the measured profiles quite well, including the twin off-center peaks and a local minimum at 

4 of 17 


American Institute of Aeronautics and Astronautics Paper 2010-4001 



the centerline, z=0. The twin off-center peaks shown in these figures indicate that the most unsteady or energetic part of 
the exit jet occurs away from the centerline. Note that the maximum levels for the perturbation velocities for the blocked 
feedback channel configuration are approximately 20-25% of the corresponding mean velocity levels, which is typical 
of free jets. 24,25 On the other hand, the perturbation velocities for the baseline 124D sweeping jet are almost twice as 
large, which indicates that the sweeping jets are more energetic in nature, and are accompanied by higher levels of flow 
unsteadiness. 

Next we show the evolution of velocity profiles with the downstream distance x for the baseline 124D case in Fig. 8, 
where the computed time-averaged and perturbation velocities are compared with the experimental data at x= 6, 12, and 
18 mm stations. As expected, the peak velocities drop in magnitude with increasing value of downstream distance x, and 
the velocity profiles become fuller by spreading out in the spanwise direction. The simulations over-predict the average 
velocities and under-predict the perturbation velocities by approximately 10-15% compared to measured data. The 
largest differences are observed at the x=18 mm location, where the experimental data displays significant asymmetry 
around centerline. Even though there are quantitative differences in computed and measured velocities, it is important to 
note that the computations predict the correct general trends, which can be very useful for conducting parametric studies. 

The final set of results for the 124D configuration examines the effect of varying the mass flow rate flowing through 
the actuator. Experimental data at x=6 mm location is also available for the mass flow rates of 0.010 and 0.020 lb/sec 
for this configuration. The numerical simulations indicate that both of these mass flow rates produce sweeping jets. The 
velocity profiles showing the effect of varying the mass flow rates are shown in Fig. 9 at x= 6 mm. As seen in earlier 
results, the computed profile shapes are fuller compared to experimental results. The average velocities are in fairly 
good agreement with the measurements except for the mass flow rate of 0.020 lb/sec. The perturbation peak velocities 
for all three mass flow rates are under-predicted by approximately 10-20% compared to measured values. However, 
the general shape of the perturbation velocity profiles and trend with variation in mass flow are similar to the measured 
data. Once again, the biggest differences are observed for the higher mass flow rate (0.020 lb/sec) case. It is not clear if 
compressibility effects or some other factors at the higher speeds encountered for the larger mass flow rate cause these 
differences. 

B. Type I, 90D actuator 

Computations were also performed for the actuator with 90° divergence angle for a mass flow rate of 0.015 lb/sec. The 
time-averaged velocity profiles are compared with the experimental data in Fig. 10 at a distance of 6 and 12 mm from 
the exit plane. The velocity peaks for this case are slightly over-predicted. Overall comparisons at the 6 mm location 
are quite good. However, at the 12 mm location, the experimental data displays an asymmetric plateau region in the 
mid-section for mean profiles, a feature that has not been observed in any of the results for the 124D configuration. The 
cause of such different behavior is not clear at this point, and it might be worth repeating the experiment for this case. 
Except in this mid-section region, the agreement between the computed and experimental results is quite good. 

The next set of results shows the effect of decreasing the divergence angle of the actuator from 124 to 90 degrees. 
Fig. 1 1 (a) shows the comparison of time-averaged velocities normalized by their respective peak levels for the 90D and 
124D configurations at a distance of 6 mm from the exit plane. The velocity profiles are broader in shape for the 90D 
case compared to the 124D case, which is indicative of larger jet width. In addition, based on the results shown in Fig. 1 1 
(b), where the peak mean velocities are plotted as a function of distance x, it is clear that a reduction in the divergence 
angle of the actuator causes the peak velocities to decrease. The computational results for the time-averaged velocities 
are in general higher for both configurations compared to the measured values. Except for the station closest to exit 
plane (x=3 mm), the measured and experimental curves are nearly parallel to each other indicating good agreement for 
decay rates of mean velocities with downstream distance x. It is conjectured that measurements at the x= 3 mm location, 
which is very close to the actuator could suffer from interference effects between the hot-wire probe and the actuator 
surface. It is encouraging to note that the effect of varying the divergence angle on peak levels is also well predicted for 
mean velocities as indicated by the incremental shift in the curves shown in Fig. 11 (b). 

C. Type II, curved-wall actuator 

The next set of results presented here correspond to the Type II, curved- wall actuator that was shown in Figs. 2-3, and 
belong to the class of devices that are currently being used in practical flow control applications to mitigate separated 
flow regions. 6,7 Time-averaged and perturbation velocities for this configuration are presented in Figs. 12(a) and (b) at 
x= 6 mm, where the computed and measured velocity profiles are compared. Results from the 124D case discussed in 
the earlier section are also added to these figures for comparison. The velocities presented in these figures have been 
normalized by the corresponding maximum time-averaged velocities for each case. Several observations can be drawn 
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from the results shown in Fig. 12. First of all, the computed profile shapes are in good agreement with their experimental 
counterparts for both the average, and the perturbation (rms) velocities. The normalized perturbation velocities are 
under-predicted by approximately 10% compared to measurements. The computations also predict the correct trends 
between the 124D and curved- wall actuators. For example, the mean and perturbation velocities for the curved- wall 
actuator are much broader in shape compared to the 124D actuator. Also note that as a fraction of mean velocities, the 
perturbation velocities for the curved- wall actuator are nearly twice in magnitude when compared to the 124D case. The 
larger spanwise footprint and higher perturbation velocities produced by the curved actuators potentially make these 
more effective devices for flow control applications. 

Note that for the Type II curved-wall actuator, the velocity profiles display different trends when compared to the 
Type I actuators that were formed with straight edges and sharp comers discussed in previous figures. For example, the 
twin off-center peaks and a sizable dip in the middle is observed in the velocity profiles for the curved-wall actuator in 
contrast to the single velocity peak at the jet centerline for Type I actuators discussed earlier in the paper. To understand 
the reason for such significant differences in jet characteristics of these actuators, the experimentally measured time 
history of the velocity magnitude at z=0 and 10 mm, and x=6 mm for the baseline Type I, 124D actuator is compared 
with the corresponding results for the Type II, curved- wall actuator in Fig. 13. The velocity variation with time for 
these two type of actuators display very different patterns. For the 124D actuator, the jet velocity at centerline (z=0) 
spends considerable time at the higher levels before dropping in magnitude. In contrast, the jet velocity at centerline 
for the curved- wall actuator stays at the higher values very briefly, then drops in magnitude and stays at low values for 
extended periods. On the other hand, at the off-center (z=10 mm) location, the velocities for the curved- wall actuator 
spend approximately equal time between the high and low values as seen in Fig. 13 (b). Also note that the peak velocity 
levels for the curved actuator at z=10 mm are almost twice in magnitude when compared to the 124D actuator results. It 
is apparent from these velocity signals that the curved- wall actuator jet spends more time at the outer edges of the active 
region and passes very quickly over the centerline, thus creating twin peaks at the two ends and a dip in the velocity near 
the centerline. 

Next, the internal and external flow fields associated with the fluidic actuators are examined to improve our under- 
standing of the underlying flow physics for such devices. For the curved- wall actuator under consideration here, the 
instantaneous velocity magnitude contours from numerical simulations are shown in Figs 14 at the mid-plane (y=0) 
passing through the computational domain. The velocity vectors are also shown in these figures to delineate the instan- 
taneous flow direction. Fig. 14 (a) corresponds to the time when the jet is aligned with left side of the exit wall. Note 
that when the jet aligns with the left side, it forces some of the fluid to enter the right side of the feedback channel and 
creates a low speed region on the left side. The curving of the jet in this position increases the pressure in the left section 
of the chamber, which pushes the jet towards the right side. The exact opposite flow pattern develops inside the chamber 
when the jet aligns with the right side at a different time during the simulation cycle as shown in Fig. 14 (b). Such a 
cyclic, bi-modal process sets up an oscillatory jet emanating from the nozzle exit, which in turn can energize the outer 
flow and enhance the aerodynamic performance of control surfaces, as demonstrated in Ref. 7. 

Finally, we examine the behavior of the peak velocities as a function of distance x, measured from the actuator exit 
plane for the curved- wall actuator. These comparisons for the time-averaged mean velocities are presented in Fig. 15 
for x varying from 3 to 24 mm. The results for the baseline 124D actuator are also shown in this figure for comparison 
purposes. The computed velocities for the curved- wall actuator are substantially smaller compared to the measured data, 
the largest differences occurring near the exit plane. The differences become smaller with increasing values of x. Note 
that the peak velocity magnitudes for the 124D case are much larger compared to the curved actuator case for the same 
mass flow rate. However, as seen in the Fig. 12, the curved-wall actuator produces a much fuller profile that spreads 
over a larger spanwise extent in the field. The CFD results presented here display the correct qualitative and quantitative 
trends in velocity profiles with variation in geometric shape for the two different type of actuators examined here. 

D. Jet characteristics 

From the flow-control perspective, it is important to examine and understand key properties of the jet flow produced by 
the AFC devices. These include the frequency and width of the oscillating jets. The frequency plays an important role 
in how quickly the jet mixes with the outer flow, whereas the jet width is representative of the extent of the physical 
domain occupied by the most energetic and unsteady part of the outflowing jet. Such quantitative parameters are helpful 
in conducting parametric studies to compare the relative efficiency of different AFC devices. Effect of geometric and 
input parameters on these quantities are discussed in the following section. 

During the simulations performed with the PowerFLOW®code, significant amount of time-dependent data were col- 
lected in pre-selected regions and at probe locations to construct the time-series, and to extract the oscillation frequency 
using spectral analysis. In the experimental work, hot-wire data for instantaneous velocities were collected along pre- 
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selected rake locations to extract statistical information. A close examination of the spectra obtained from the unsteady 
data indicates the existence of a primary frequency that contains the most significant energy content, and hence it is the 
most important frequency for comparing relative characteristics of such devices. The primary frequency of oscillations 
at x=6 and z=10 mm probe location obtained from the experimental data and simulations are compared in Fig. 16 for 
the 124D and curved-wall actuators corresponding to mass flow rates of 0.010, 0.015 and 0.020 lb/sec. As expected, 
the oscillation frequency increases with the mass flow rate due to increase in flow velocities. The predicted frequencies 
are somewhat higher compared to their measured counterparts, but replicate the experimentally observed trends with 
the variation in mass flow rate fairly well. Except for the lowest mass flow rate (mdot=0.010 lb/sec), the simulated and 
measured frequencies for the curved- wall actuator are slightly lower compared to the baseline 124D case. 

Final results presented here compare the growth of jet half-width with the downstream distance x. The jet half-width 
is defined as the distance between two points where the mean velocity drops from the maximum to half of that value. 
For the curved- wall actuator that produce twin peaks in the velocity profile, we take the average of the left and right side 
jets. The jet half-widths for the 124D and curved-wall actuators are plotted as a function of distance x in Fig. 17. The 
CFD results for the jet width are in excellent agreement with the measured data. The computational results also correctly 
predict the increment in jet width due to configuration change from the 124D to the curved- wall actuator. 

VI. Concluding Remarks 

Numerical simulations have been performed using the PowerFLOW® code to compute the unsteady flow field 
generated by a family of fluidic actuators. These simulations provide useful information about the flow physics in 
the internal and external regions of the fluidic actuators. The bi-stable feedback mechanism that produces sweeping, 
oscillatory jets is demonstrated with the help of computed flow field plots. The present results also demonstrated that 
blocking the feedback channel inhibited the sweeping action of outflowing jets. 

The effect of varying the inflow mass flow rate, as well as key geometric shape parameters on the jet velocity profile 
and spanwise spreading rate were presented. Overall agreement of the computed time-averaged mean and perturbation 
velocities with measurements for the Type I, sharp-edged actuators was quite good. For the curved-wall actuators of 
Type II, the mean and perturbation velocities were under-predicted by about 25% near the exit plane, the agreement 
improving with downstream distance. The simulations replicated the correct overall trends and profile shapes when 
compared with the experimental data. 

The primary frequency of oscillation of the sweeping jets produced by these fluidic actuators were extracted from the 
unsteady data using spectral analysis. The computed frequencies compare favorably with measured values and displayed 
the correct trends with respective variation of input mass flow rate. Results were also presented for jet half-width for the 
baseline mass flow rate of 0.0150 lb/sec, where the computed jet widths for both the sharp-edged 124D actuator and the 
curved- wall actuator were found to be in excellent agreement with the measured data. It was observed from these results 
that for the same inflow mass flow rate, the Type II fluidic actuator with curved- wall interior geometry produced much 
higher spreading rate for the outflowing jet, and is potentially a better performing device for flow control applications. 

Future work will focus on simulating the effect of a series of fluidic actuators on external flow over an aerodynamic 
configuration of practical interest to examine the effectiveness of such AFC devices for separation control. Parametric 
studies will be conducted to determine the optimum spacing between actuator arrays to achieve the desired control 
authority with minimum mass flow. An effort will be made to investigate the feasibility of simultaneously simulating the 
installed actuators or formulating suitable boundary conditions to approximate the effect of such devices on the external 
flow field in order to improve computational efficiency. 
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Figure 1. Schematic of Type 1, 124D fluidic actuator, top view. 




Figure 2. Schematic of Type II, curved actuator, top view. 


9 of 17 


American Institute of Aeronautics and Astronautics Paper 2010-4001 



Figure 3. Experimental set-up for Type II curved actuator. 



Figure 4. Schematic of PowerFLOW® computational domain for Type I fluidic actuator. 
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Figure 5. Total velocity contours for Type 1, 124D actuator, blocked feedback channel. 
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Figure 6. Total velocity contours for Type 1, 124D actuator, original configuration. 
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(a) time-averaged velocities (b) perturbation velocities 

Figure 7. Effect of blocked feedback channel on velocity profiles for Type 1, 124D actuator, x=6 mm. 




(a) time-averaged velocities (b) perturbation velocities 

Figure 8. Evolution of time-averaged and perturbation velocities for Type 1, 124D actuator. 
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(a) time-averaged velocities (b) perturbation velocities 

Figure 9. Effect of varying mass flow rate on velocity profiles for Type 1, 124D actuator, x=6 mm. 



Figure 10. Time-averaged velocity comparisons for Type I, 90D actuator. 
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(a) mean velocity profiles (b) decay of mean velocities with x 

Figure 11. Effect of divergence angle on mean velocities for Type I actuators, x=6 mm. 




(a) time-averaged velocities (b) perturbation velocities 

Figure 12. Normalized velocity comparisons for Type 1 124D and Type II curved-wall actuators at x=6 mm, mdot=0.015 lb/sec. 
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(b) z = 10 mm 

Figure 13. Time variation of velocities at x=6 mm, experimental data. 
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Figure 14. Total velocity contours for Type II, curved-wall actuator, m=0.015 lb/sec. 



Figure 15. Decay of peak velocities with x-distance for Type II, curved actuator, mdot=0.015 lb/sec. 


16 of 17 


American Institute of Aeronautics and Astronautics Paper 2010-4001 


400 - 


350 - 


N 

X 

> 

£ 300 

CD 
O" 

<D 


250 - 


'/ 


/ 

✓ / 
x /' 

/ ,/' 

X 


A 

/ 


/ 


// 


/ 


// 

'/ 


/ 


/ 


~V~ 


EXP (124D) 
EXP (curved) 
CFD (124D) 
CFD (curved) 


200 - 


_L 


_L 


0.01 


0.015 

Mass flow rate (Ib/sec) 


0.02 


Figure 16. Frequency variation with mass flow rate for 124D and curved-wall actuators. 



Figure 17. Comparison of jet half-widths for 124D and curved-wall actuators, mdot=0.015 lb/sec. 
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